Excitonic effects in the optical conductivity of gated graphene 
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We study the effect of electron-electron interactions in the optical conductivity of graphene under applied gate 
and derive a generalization of Elliott's formula, commonly used for semiconductors, for the optical intensity. We 
show that excitonic resonances are responsible for several features of the experimentally measured mid-infrared 
response of graphene such as the increase of the conductivity beyond the "universal" value above the Fermi 
blocked regime, the broadening of the absorption at the threshold, and the decrease of the optical conductivity 
at higher frequencies. 
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Since its isolation in 2004 |Q], the demonstration of the 
presence of electron-electron interactions in the electronic 
properties of graphene has been elusive [2 |. The first clear 
manifestation of interaction effects comes from the recent 
measurement of the fractional quantum Hall effect at filling 
1/3 EIIU. To understand the impact of many -body physics in 
graphene and disentangle it from purely band structure effects 
it is essential to fully understand graphene's electronic and 
transport properties. Furthermore, because Coulomb interac- 
tions can be controlled externally by choosing the dielectric 
constant of the substrate on which graphene is deposited, it 
is possible to tune the electron-electron interactions in a way 
that cannot be done in most materials. 

An obvious property where to look for correlation effects 
is the optical response. It is well established that in ordinary 
semiconductors electron-electron interactions are essential to 
explain their optical properties [5|. For energies below the 
gap, a particle and a hole created by the excitation process 
can form a bound state (an exciton), which appears as a well 
defined peak in the absorption spectra, in the region within 
the energy gap. For energies above the gap, a renormalization 
of the band and an enhancement of the optical conductivity is 
observed and attributed to the Coulomb-mediated scattering 
between the electron-hole pair. These effects are embodied 
in the famous Elliott's formula [5| that describes the intensity 
of optical absorption close to the gap edge. This is true for 
traditional semiconductors in both three and two dimensions 
@. As we are going to show the situation in graphene is, 
however, quite different. 

Firstly, in graphene there is no gap between the lower 
(it) and upper (w*) bands, the effective theory describ- 
ing graphene's low energy physics is the massless two- 
dimensional Dirac equation, and the Coulomb problem in this 
case has no true bound states, but resonances flTl flOl . These 
facts prevent the formation of bound excitons in graphene. 
On the other hand, exciton resonances can exist and a signa- 
ture of their presence should be seen in its optical properties 
lfTTT4l31 . Although the system is gapless, the possibility of 
doping it by using a gate creates an energy interval, of width 
2/i (/i is chemical potential), where light absorption is forbid- 
den due to Pauli's exclusion principle. In tradicional semi- 



conductors, such a regime is known as the Moss-Burshtein 
effect. All theoretical investigations of the optical conductiv- 
ity of gated graphene developed so far evaded the problem of 
electron-electron interactions and the discussion of excitonic 
effects [ 1 6—20 1 . However, a detailed analysis of the experi- 
mental data lfT2lfT4l shows that graphene's optical conductiv- 
ity deviates considerably from the non-interacting prediction 
H21Q3L e yen with the inclusion of band structure corrections 
to the Dirac spectrum ETI . disorder, and phonons |fl8l . 

The alluded deviations occur in both neutral and doped 
graphene, but they are far more evident in the latter case. For 
neutral graphene, in the visible range of the spectrum at ener- 
gies around ~ 3 eV, the measured percentage of transmitted 
light, T, is smaller than the value T = 1 - ira ~ 97.7% (a 
is the fine structure constant), predicted by the non-interacting 
theory [12|. Band structure effects alone cannot account for 
the measured deviations ll2Tl . In the left inset of Fig. [TJit is 
shown that calculations based on ab-initio methods can de- 
scribe, qualitatively, the data (see also Ref. [22b . 

For gated graphene, experiments [ 14] show even more no- 
ticeable deviations from the non-interacting picture, which 
predicts that the optical conductivity of graphene has the form: 
cr(bS) — o-Q9(fku — 2/i), where ctq — ire 2 /(2h) is the so-called 
"universal" AC-conductivity, hui is the photon energy, and 
9(x) is the Heaviside step function. The deviations seen in the 
data are of five different types: (f) finite absorption below 2/i, 
which is due to both inter-band and intra-band elastic and in- 
elastic scattering processes; (if) broadening of the absorption 
edge around the energy threshold 2/i; (fff) an enhancement of 
the conductivity above the universal value, <7o, in the energy 
range between 2/i and 2/i + E* , where E* is a characteristic 
energy scale; (fv) a reduction of the conductivity bellow o-q, at 
energies above E* , with the conductivity, as a function of fre- 
quency, having a positive curvature; (v) the imaginary part of 
the conductivity is larger than the value predicted by the non- 
interacting model for energies hut 2/i. Finally, we point out 
that the optical conductivity curves, for diferent gate voltages 
O, collapse on top of each other, when re-plotted as func- 
tion of oj/ [i, implying that the mechanism causing deviations 
from the non-interacting approximation must be intrinsic. The 
analytical microscopic theory we develop in this letter, which 



2 



includes excitonic effects, accounts for items (if), (Hi), (iv), 
and (v), and it also partially accounts for item (i), although not 
completely, since we have not included intra-band scattering 
in the calculations. 

The central result of this work is the calculation of the con- 
ductivity of doped graphene, that is, at finite gate voltage, V g . 
The calculation takes into account electron-electron interac- 
tions at the exchange level. Its important to stress that the 
theory given below has no adjustable parameter (except for 
the use of an effective temperature -see discussion below); 
the density of impurities is determined by fitting the DC con- 
ductivity, that is, it is fixed a priori to the calculations of the 
optical conductivity. The final result for the latter quantity 
reads: 

(1) 

where all quantities, to be given below, depend explicitly on 
oj, and on temperature. In Eq. ([T]), e > is the electron 
charge, Vp is the Fermi velocity, d x xx , fc is the matrix element 
of the dipole operator with A, A' = c, v referring to the it* 
(it) band index, respectively, x'xx' k * s trie optical suscepti- 
bility of graphene in the independent electron approximation, 
and r>A'.fc is the vertex correction due electron-electron in- 
teractions. If we assume r^v.fc = 1, Eq. (|T|) gives the 
"universal" value for the conductivity of graphene, Go- The 
form of the vertex function reads: T\y ^ = [1 + St )Cl fc 1 i] _1 > 
where s VCf k,i is defined below. This result shows that if 
rAA'.fc < 1, its contribution to Eq. ([T]) acts to reduce a rel- 
atively to (To. On the other hand, when T\x',k > 1 it con- 
tributes to the enhancement of a over <tq. We show below 
that the effect of the vertex correction on a depends on lu. 
The conclusion is then: the importance of electron-electron 
interactions in graphene can manifest itself in the behavior 
of a by studying how the data deviates from ctq. More- 
over, the optical conductivity relates to the transmittance by: 
T = [(1 + 5Rcr/(2e c)) 2 + (3a-/(2e c)) 2 ]-\ where e is the 
vacuum permittivity and c is the speed of light. Our results 
are summarized in Fig. [T] In the left panel, we depict the real 
part of the conductivity: the dashed black curve is the non- 
interacting prediction, considering the effect of disorder; the 
solid red curve is the prediction taking into account excitonic 
effects and the same amount of disorder. The right panel refers 
to the imaginary part of the conductivity. In both panels the 
experimental data [14| is also shown. 

The five items mentioned previously relatively to the data, 
can easily be observed to be present in the theoretical curve; 
(i): For Huj < 2/x, there is a finite conductivity, where the 
non-disordered independent electron model predicts zero ab- 
sorption. This modification over the naive result is mainly a 
consequence of disorder, and excitonic effects play no role, 
as it should be. The discrepancy of about a factor of two 
between the theory and the data, is to be expected and is a 
strength of the model, since additional intra-band scattering 
(the model only includes inter-band scattering) will transfer 
spectral weight from the Drude peak to this frequency region, 




(0/ p, (£>/ [l 

FIG. 1: (Color online) Real and imaginary parts of the conductiv- 
ity. The black dashed line is obtained from the independent electron 
model with disorder. The solid red line is the model given in the text, 
which includes excitonic corrections and disorder. The experimental 
curves, refer to a gate voltage of 28 V, to which corresponds a Fermi 
energy of fj, ~ 0.18 eV. In the experiment! 14 1, the chemical potential 
was given by fi = 0.03474\/V g eV. The calculation took an effec- 
tive temperature of 120 K (see the discussion in point (ii) in the text). 
In the left inset we give the transmission, T, of neutral graphene, as 
function of ftw, in comparison with the experimental data of Ref . [121 
In the right inset we show the DC conductivity as a function of gate 
voltage, comparing both the experimental data (blue line) of Ref. 
Q4]and the thoretical fit; the fit fixes the concentration of impurities, 
to the values of m = 2 x 10 11 cm~ 2 , for resonant scatterers, and 
rii — 1 x 10 11 cm -2 , for charged ones. The determined impurity 
concentrations are then used in the calculation of the optical conduc- 
tivity. Therefore, the calculation of the optical conductivity has no 
fitting parameter. The labels (i), (ii), (iii), (iv), and (v), refer to the 
five items discussed in the text. 



enhancing the optical response, (ii): It is clear that disorder 
and temperature effects combined, account for the broadening 
of the Fermi step around 2/i. As discussed in the context of the 
optical response of the biased graphene bilayer 11231 . charge 
density fluctuations, associated with the random electrostatic 
potential present at the graphene-Si02 interface, smear the 
Fermi energy, leading to and effective higher temperature than 
that measured by the finger temperature of the cryostat (in the 
experimente [ 14], the temperature of the finger was 54 K); as it 
is well known, an increase of the temperature broadens Fermi 
edge ETl . as seen in the experiments (we have used the same 
effective temperature as that used in bilayer graphene, 120 K; 
this makes sense because the random electrostatic potential al- 
luded above is a property of the substrate), (iii): For hut > 2/i, 
an enhancement of the conductivity above <7o is seen. This 
behavior extends to the characteristic energy scale E* . The 
energy interval over which Sftcr > uq is about the same in the 
theory and in the data. Finally, we should stress that the results 
of Fig. [T]were obtained from an approximate analytical solu- 
tion to the vertex functions. Corrections to this approximation 
may also account for the small differences between the data 
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and the theoretical curves, (zv): For Huj > E* , the conductiv- 
ity goes below <tq, over a large energy range, presenting a pos- 
itive curvature. Exactly the same behavior is seen in our calcu- 
lations, with the decrease of the conductivity below er having 
about the same numerical value, (v): The imaginary part of 
the conductivity, 3cr, is larger, for energies above ~ 3/i, than 
that predicted by the non-interacting model with disorder and 
fits well the data. In conclusion, Eq. ([T} does explain all fea- 
tures, (ii)-(iv), of the experimental data J14l . making clear that 
electron-electron interaction effects are manifestly present in 
the optical conductivity of gated graphene for Hu > 2/Lt. Fi- 
nally, we must note in passing that we have chosen to fit a set 
of data where no renormalization of the Fermi velocity was 
measured lfl4l . since our model does not include the renor- 
malization of the electronic velocity For all data sets such that 
V g > 28 V, no Fermi velocity renormalization was measured; 
we have chosen the first of this class of sets, namely that for 
which V g = 28 V. 

Before discussing the calculation of the electron-electron 
corrections to the optical conductivity, it is instructive to 
show that the formalism based on the polarization concept 
gives the well known result for the "universal" conductivity 
of graphene. As it will be clear, this formalism bypasses the 
calculation of the Kubo formula lfl8l \l9l . The low-energy 
Hamiltonian of graphene, with a vector potential, A(r, t), 
reads: H = vpcr ■ [p + eA(r, t)] = H + Hj . Assuming 
an electric field of the form E(t) = Eoxe~ luJt , the interaction 
of the electrons with light is given by Hj = ev F Ao<j x e~ lut = 
ev F a x Ao(t), with Aq = —iEq/lo. The polarization along the 
x— direction is: 



P x (t) = J2 / dr{i>i(r,t)dMr,t)) 



(2) 



result for the thermal average of the operator P vc ,k- 

(Pvc,h) = v F eA d x cv ^ k y :' , (5) 

where np(x) is the Fermi function, E c / vk — ±v F hk, and 
r(w) is the level broadening (a similar equation follows for 
(Pcv,k))- We note that the broadening T(ui) is not a con- 
stant, but follows a well prescribed model, where the scat- 
tering is dominated by resonant and charged scatterers I9l l24l : 



F(w) = n^ s \ F ha T '' > (u>)+n { i w, v F ha^ n> (uj), where 

is the concentration of impurities and Op S ^ ch \ui) is the ex- 
act transport cross section, for each type of scatterers. Once 
(P\\',k) is known, the calculation of the polarization per unit 
area, p x = P x (uj)/A, follows from Eq. (|4ji as: 



(«*). 



(rs/ch) 



2e 2 v F E Q ^ 

Px = ~ l — — — 2^ s ( k > 

k,X^\' 
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n F {Ex,k) - n F (E X ',k) 

uh - E\ t k + Ey,k + *r ' 

(6) 

where s(k) = sin 2 [6>(fc)], and the optical conductivity is sim- 
ply given by a = p x /Eq. Equation (|6ji for the conductiv- 
ity of graphene is the same one obtained by computing the 
Kubo formula from the current-current correlation function 
CTll . If we consider, as a simple example, the limiting case 
where T — »• 0, fj, = 0, and zero temperature, Eq. ([6]) gives 
Px/E = ne 2 /(2h) = a Q . 

We can extend the formalism above to the case where 
electron-electron interactions are considered, and use it to de- 
rive the central result of this work, Eq. ([lj. The formalism 
we present below, gives results equivalent to the solution of 
the Bethe-Salpeter equation, but following a more straight- 
forward path. Electron-electron interactions in graphene are 
described by the Hamiltonian: 



where d x is the dipole operator, d x = —v F ea x , and Vv(r, t) 
is the field operator: 



ff e _ e = - / drdr'^(r)^{r')V(r-r')i/j(r')il;(r) , (7) 



Vv(r, t) = ax,k,a(t)ipx,k( r ) > 

k,\ 



(3) 



with tpx.k{ r ) the eigenfunctions of H (2), ax,k,o the destruc- 
tion operator of an electron in band A, with momentum k = 
(k x , k y ), and spin projection a (A = 1 for A = c and A = — 1 
for A = v). The matrix element of the dipole operator reads: 
d x x y k = {\'e lB ^ +Xe- w ^)/2, withtanfl(fe) = k y /k x . 
Then, the final expression for the polarization (per unit time) 
P x (t) is, in second quantized form, given by: 



P x (t) = -Ev F Y, X>U,a(*K',fc.*(*M,V 



k ■ 



(4) 



,A,A' 



and Hj reads: 

Hi = v F eA(t) E ct ,a,A' E fc "a t -r"-\' a- -'/.\.\-.a- • 

We introduce the dipolar operator P VCyk = a} y k a c , k , and 
seek the solution of its equation of motion, —ihdP VCtk /dt = 
[H,P vc , k ]. The procedure is simple, and gives the following 



where V(r) is the Coulomb potential. Since we are studying 
the conductivity of gated graphene, the Coulomb potential is 
screened, and at the level of the Thomas-Fermi approxima- 
tion (corresponding to the static RPA limit) the Fourier trans- 
form of the Coulomb potential reads: <P2d(q) = e 2 [2eo(q + 
Qtf)}^ 1 , with the Thomas-Fermi momentum, qtf, defined 
as (7tf = 4 
area and a g = 



/im e , where n e is the electron density per unit 
2 /(Aireov F K) (for graphene on top of SiC»2 



we replace e by e esio 2 an d e 2 /(4e esio 2 t '-F 7r ^) ~ 1/2; in 
our calculations we used a g = 0.4). Introducing the field op- 
erators ([3]l in Eq. |7}, the electron-electron interaction in the 
momentum representation reads: 



q} X 
\-q,(7 \2,k' — q,a' 



E X^dC^L/m 

{Aj } k,k' ,q a, a' 

ax 3 ,k',^a\ A ,k,cxf{k' , -q)x 2 ,\ 3 f( k > > ( 8 ) 



where f{k,q) Xa ,x b = (1 + \ a \ b e- mk+q) - e{k)] ) . As be- 
fore, we seek the equation of motion for the operator P vc , k = 
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ol k a c.k, but now in the presence of electron-electron inter- 
actions. To that end, we need to compute the commutator 
of P vc ,k with H e ^ e , in addition to the one we have already 
determined for the non-interacting theory. The evaluation of 
{[H e - e , P V c,k]}, after making the usual Hartree-Fock decou- 
pling the quartic terms, gives: 



([He-e, Pvc,k\) — @ 



cxch. 



-c 



dens. 



(9) 



where the several terms have the following physical inter- 
pretation: (7 cxch - i s the exchange correction to the bare 
bands; C cxclt - are the excitonic contributions, coming from 
the Coulomb interaction among the electron-hole pair formed 
by the absorption of light; C dcns ' is the interaction of the elec- 
trons and holes in the bands with the electron gas formed by 
the gating; and C n_1 are non-linear terms in the field intensity, 
Eq, due to the interaction of the average of the operator P VC:k 
with itself and with the average of the operator P CVtk (see also 
Ref.|25). 

Let us consider the calculation of the vertex corrections, 
F\y. k , entering in Eq. (|TJ. We ignore the non-linear correc- 
tions included in the commutator {9}. The term C dcnslty can 
be shown to give a zero contribution to the polarization of the 
material. From the commutator the full equation of mo- 
tion for the operator P VCt k gives: 

(-wh - E vh + E c _ k )(P VC:k ) = E -^dc V k Sn Fik 

ltd 

+ tin Fth -^-T^2<P2T>(q)(Pvc,k-q)F + , + (k,q) 



q.a 



+ 5n F , k ^^2<p 2 D{q)(Pcv,k-q)F--(k,q), (10) 



with 5n F/h = [n F (E C;h ) ~ n F (E v>k )] and F A)A (fc,g) = 
1 + Acos[#(fc) — 9{k — q)]. Eq . ([To| and the equivalent one 
for (P cv k)< obtained from Eq. ( 10 1 by interchanging the la- 



bels c and v, form a set of two coupled integral equations. The 
approximate solution of this set of equation is obtained as fol- 
lows: since we are interested in the linear response, we write 
Pvc,k = EoXvc.k- Additionally, the function Xvc,k is written 
in terms of the x®c k (which represents Xvc.k in the absence of 
electron-electron interactions) as Xvc,k — ^vc,kXvc fc' where 
r„ c ,fe is the vertex function. After these operations, a set of 
coupled integral equations for the vertex functions is obtained 



from Eq. ( 10 1, which is then solved by the method of Pade ap- 
proximantslO . For simplicity, we have only considered the 
leading approximant. The method of solution described here 
is equivalent to summing an infinite sub-set of Feynman dia- 
grams, and it gives the analytical solution for the vertex in the 
form T A A',fc ~ [1 + s„ Cifcj i] _1 with the s AA ',fe,i defined as 



where, 
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UJ 



v F esmu k 



©VA.fc- (ID 



e 



± 

XX' k' 



1 
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This procedure for determining the vertex function is equiva- 
lent to the RPA solution of the Bethe-Salpeter equation [26 1 ; 
we believe, however, that our method is more straightforward. 
Equations (Hi and ( 12 1 can be used to determine the vertex 



d6 q I qdqip 2 v{q-k)x X y, q [^±cos{0 k -6 q )\. 

(12) 



function in closed form. The bare susceptibility, x^w k> I s 
given by 



X°xy 



sin Uk - 



ev F fn F (Ex,k)-n F (Ex' 



\hw~E. 



X'k 



-E 



k) +x ,0{k) 



X,k 



-iT 2v F hk 



(13) 



where 9(k) is the Heaviside step function. Unfortunately the 
integrals cannot be evaluated analytically but can easily be 
obtained by numerical integration. 

In summary, we have developed a general formalism for 
taking into account excitonic effects in graphene. These ef- 
fects were computed taking the lowest Pade approximant, but 
higher order ones can be easily introduced and may improove 
the fit of the data. Additional many-body effects, such as in- 
elastic phonon scattering can also be easily included, which 
makes the formalism a powerful one to study the combined 
effects of disorder, electronic correlations, and phonon scat- 
tering in graphene's optical spectrum. Extending the formal- 
ism to visible range of the spectrum is also feasible, by com- 
puting the eigenproblem of the free theory including trigonal 
warping effects. Finally, we have showed that our formalism 
solves discrepancies between the data [ 14 1 and the theoretical 
results provided by the independent electron picture. 
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